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Abstract 

In this paper we study the cross section at leading order in 1/Q for polarized 
Drell-Yan scattering at measured lepton-pair transverse momentum Qt- We 
find that for a hadron with spin 1/2 the quark content at leading order is 
described by six distribution functions for each flavor, which depend on both 
the lightcone momentum fraction x, and the quark transverse momentum fc^. 
These functions are illustrated for a free-quark ensemble. The cross sections 
for both longitudinal and transverse polarizations are expressed in terms of 
convolution integrals over the distribution functions. 
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I. INTRODUCTION 



The measurements of unpolarized structure functions in deep inelastic scattering (DIS) 
of leptons off nucleons and nuclei and those of polarized structure functions in scattering 
of longitudinally polarized electrons off longitudinally polarized nucleons |l| have yielded 
the lightcone momentum distributions fx(x) for quarks in various targets and the helicity 
distributions g\(x) in protons and neutronsQ. These measurements, and particularly their 
interpretation, have shown the importance of understanding the relation of these distri- 
butions to the structure of the target. The distributions fi(x) and g\(x) characterize the 
response of the hadron in inclusive DIS at leading order in the transferred momentum Q. 
In inclusive deep inelastic lepton-hadron (£H) scattering the quark transverse momentum 
is not observable, since it is integrated over. In the Drell-Yan (DY) process at measured 
lepton-pair transverse momentum Qt, however, quark transverse momentum does enter in 
observables, notably in the angular distribution of the lepton pairs. The main point of this 
paper is the discussion of quark transverse momentum in polarized Drell-Yan scattering. We 
will restrict ourselves to leading order and discard contributions which are suppressed by 
orders of 1/Q. We will also not discuss QCD radiative corrections, giving rise to logarithmic 
corrections. 

For inclusive deep inelastic £H scattering, assuming only one flavor, the hadron tensor 
is given as the imaginary part of the forward virtual Compton amplitude, for large virtual 
photon momentum q (Q 2 = —q 2 large) given by the sum of the quark and antiquark handbag 
diagrams of Fig. [L| The basic object, encoding the soft physics of the quarks inside the 
hadron, is the correlation function [^,[| 

^(PS;k) =/t04 e^(PS\^(0)Ux)\PS)c, (1-1) 

where k is the momentum of the quark. The vectors P and S are the momentum and 
spin vector of the target hadron. Evaluating the hard part, the scattering of the virtual 
photon off the quarks, it turns out that the structure functions in the cross section become 
proportional to /i(xbj) and gi(x^), where x b j = Q 2 /2P ■ q. The function fx is given by 

h{x) = l -J dk-d 2 k T Tr [ 7 +$ (PS; k)] , (1.2) 

where x = k + /P + . It can be interpreted as the longitudinal (lightcone) momentum distri- 
bution of quarks. The function gi appears as 

\g x {x) = ~J dk~d 2 k T Tr [ 7 + 75$0PS; k)} , (1.3) 

and can be interpreted as the quark helicity distribution in a longitudinally polarized nucleon 
(helicity A = 1). The functions fi and g\ are specific projections of $. Which projections of 



1 Another often used notation is q(x) for the lightcone momentum distribution and Aq(x) for the 
helicity distribution (q = u,d,s, . . .). 
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$ contribute in hard scattering processes in leading order can be investigated by looking at 
the operator structure, including the Dirac and Lorentz structure, of the correlation function. 
Such an analysis requires some physical constraints on the range of quark momenta. The 
analysis of $ (integrated over k~ and h?) shows that there is one more leading function, the 
transverse polarization or transversity^ distribution h\. It is related to a bilocal quark-quark 
matrix element through || 

hi(x) = -J dk'd 2 k T Tr [ia i+ ^(PS] k)] (i = 1, 2), (1.4) 

which shows that hi can be interpreted as the quark transversity distribution in a trans- 
versely polarized nucleon. This is a chiral-odd distribution, which is not observable in 
inclusive £H scattering. It needs to be combined with some other chiral-odd structure, 
e.g., the fragmentation part in semi-inclusive leptoproduction of hadrons or the antiquark 
distribution part of DY scattering . 

In this paper we discuss one possible way to extract more information from the cor- 
relation function $. We are after the dependence on the transverse momentum fc^. One 
way to study this dependence is the observation of a hadron in the outgoing quark jet, e.g., 
in semi-inclusive £H scattering |J. This process, however, also requires consideration of 
the fragmentation functions. In this paper we study the process that is sensitive to intrin- 
sic transverse momentum and involves only quark distribution functions, namely, massive 
dilepton production or the Drell-Yan (DY) process ]10| . 



About fifteen years ago Ralston and Soper (RS) published a pioneering paper || on the 
polarized Drell-Yan process. Because we take it as our starting point, we briefly sketch its 
content. RS write down a covariant expansion for / dk~§, which is the quantity that is 
relevant in the hadron tensor for the DY process. To determine this expansion they use 
symmetry arguments and an infinite-momentum-frame analysis. They find five independent 
distribution functions, divided in one momentum probability distribution V(x,k^), two 
functions describing the quark helicity, and two describing its transverse polarization. With 
these they calculate the polarized Drell-Yan cross section with the virtual photon transverse 
momentum Qj< = yjcfe put to zero. In that case they are sensitive to four of the five 
distribution functions. When they integrate over the transverse momentum they are only 
sensitive to three distribution functions. 

We extend on these results in two ways. First, we show that RS left out one transverse 
momentum distribution, needed to describe the quark transverse polarization. This addi- 
tional function is obtained using general symmetry arguments. It also shows up in a model 
that we are going to employ later and that describes a gas of free partons. Our second 
extension is the calculation of the polarized DY cross section without constraints on Qr 
[other than it being of 0(A)], thereby becoming sensitive to all six distribution functions. 



2 The authors of Q use the name 'transversity distribution' in order to make clear that a quark of 
definite transversity is not in an eigenstate of the transverse spin operator but of the Pauli-Lubanski 
operator projected along a transverse direction. The authors of object to this nomenclature, 
because of the pre-existence of the term, and prefer to call it transverse polarization distribution. 
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We end this introduction with a remark on possible QCD corrections affecting transverse 
momenta and factorization. A difficulty of the extra scale Qt is the Sudakov effect. Soft 
gluon radiation gives rise to radiative transverse momentum. However, the large logarithms 



connected with this effect can be summed and exponentiated to Sudakov form factors |TT 
From these it becomes clear that if Qt is sufficiently low, i.e., of hadronic scale A, as com- 
pared to Q, the transverse momentum governing the process is predominantly intrinsic. 
Factorization means that the process can be written as a convolution of renormalized distri- 
bution functions and a perturbatively calculable short- distance part. For polarized DY at 



measured Qt A factorization has not been proven yet [0,0. We will not further address 



this problem here, but use the diagrammatic expansion proposed by Ellis, Furmahski, and 
Petronzio (EFP) @ to study the DY process. In this diagrammatic expansion Green func- 



tions appear, incorporating the long-range QCD physics. These correlation functions are 
connected by ordinary Feynman graphs with quarks and gluons, the hard scattering piece. 

The outline of this paper is as follows. In Sec. |IJ we give the one-photon exchange 
picture for massive dilepton production. We specify the notation in a frame where the 
two hadrons are collinear, and the axes are given with respect to which the lepton angles 



are defined. In Sec. |TTT| we analyze the quark correlation function, and find six leading 
distribution functions. In Sec. [TV] we discuss the free-quark ensemble as an example. In 
Sec. |V| we calculate the leading-order hadron tensor and present cross sections for various 
combinations of polarizations. We end with a discussion of these results. 



II. THE DRELL-YAN PROCESS 

In this section we want to discuss the cross section, kinematic aspects, and structure 
functions, for polarized Drell-Yan scattering. For a complete overview we refer for the 



unpolarized process to Lam and Tung [15], and for polarized DY scattering to Donohue and 



Gottlieb [nj, who make use of the Jackob-Wick helicity formalism. 

A. The DY cross section 

We consider the process A + B — > £ + £ + X , where two spin-^ hadrons with momenta 
P% and Pg interact and two outgoing leptons are measured with momenta k± and k%- The 
leptons are assumed to originate from a high- mass photon with momentum q = k\+k2, with 
Q 2 = q 2 > 0. We consider the case of pure incoming spin states, characterized by the spin 
vectors S 1 ^ and i.e., S\ = S% = —1. In the deep inelastic limit Q 2 and s = (Pa + Pb) 2 
become large compared to the characteristic hadronic scale of order A 2 ~ 0.1 GeV 2 , while 
their ratio r = Q 2 /s is fixed. The phase space element for the lepton pair can be written as 
d 4 qdQ, where the angles are those of the lepton axis in the dilepton rest frame with respect 
to a suitably chosen Cartesian set of axes. The cross section can be written as 



d 4 qdn 2s Q 4 ' 

where the lepton tensor is given by (neglecting the lepton masses) 

= 2 kf{k v 2 + 2 k^k\ — Q 2 g^ u , (2.2) 
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and the hadron tensor can be written as 

d A x 



W^(q;P A S A ;P B S B ) 



(2u) 



e^ x {P A S A - P B S B \ [J"(0), r{x)\ \P A S A - P B S B ). (2.3) 



Since the lepton tensor ( |2.2j ) is symmetric in its indices, we will from now on only consider 
the symmetric part of W^. 



B. Kinematics 

We define the transverse momentum of the produced lepton pair in a frame where the 
hadrons are collinear, with the third axis chosen along the direction of hadron A. One 
has q\ = Q\ < A 2 . It is convenient to work in a lightcone component representation, 
p = [p~~,P + ,Pt] with p^ 1 = (p° ±p 3 )/\/2- The momenta of the hadrons and the virtual 
photon in a collinear frame take the form 



Pa 
Pb 

q = 



2P+' A '° T 



x A Ml kQ 



,0 t 



P, 



B ■ 



2P, 



,0i 



B 



x B P B ,x A Pj[, qr 



_V2kQ' V2x a 

Q kxbMI 
V2kx b V2Q 
Q kQ 



On 



s/2k V2 



(2.4) 
(2.5) 
(2.6) 



neglecting corrections of order 1/Q 2 , indicated here and further on by an approximate equal. 
The parameter k fixes the collinear frame. One has k = x A M A /Q for the frame in which 

hadron A is at rest, k ■■ 



ix A /x B for the hadron center-of-mass frame, and k = Q/x B M B 
for the frame in which hadron B is at rest. The following Lorentz- invariant relations hold: 

' 1 = TTr ~ — ~ T— — (2-7) 

x B = w -^-^ w - " ^ , (2.8) 

(2.9) 



q + 


Q 2 


Pb 


■ q 


" PI ~ 


2P A -q 


Pb- 


Pa 


_ V~ ~ 


Q 2 


Pa 


■ q 


Pb ~ 


2P B -q 


Pa- 


Pb 


i 2P+P B 


Q 2 

X A X B 







The above relations also show that all dot products for any pair from the vectors q, P A , and 
P B , are of order Q 2 . As compared to this, the hadron momenta are almost lightlike. We 
can define the exactly lightlike vectors that in a given collinear frame have the form 



n + = [0,k,O t ], 
n_ = [k-\0,O t ], 

satisfying n + ■ n_ = 1. Given an arbitrary four- vector a, and the projector 



(2.10) 



LIU liu it V V U 

gtj> = g — n+n_ — n + n_, 



(2.11) 
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we define the spacelike transverse four-vector = 'civ, or, in coordinates in a collinear 
frame, = [0, 0, ay]. Note that for any transverse vector one has 



Q-t ■ Pa — o>t ■ Pb 



0. 



For the analysis of the hadronic tensor which satisfies g M W^ v = q u = 
tant to construct vectors that are orthogonal to q. We will use the projector 
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hi/ 



q H q u 



(2.12) 
0, it is imp or- 

(2.13) 



for this, and define 



q A 



(2.14) 



As q is timelike in DY scattering, it is useful to define a set of Cartesian axes. The in- 
direction, known as the Collins-Soper axis |T7|], is chosen as in 0, but our X- and F-direction 
are opposite. To be precise, we use (e 0123 = 1) 



Pb-Pa 
Pb-Z ^ 



Pb-Pa 
1 



P A -q 

Pa-Pb 
Pa-Z = 



P B -q 
Pb-Pa 



Pl- 



P A -q 



Pa-P 



B 



B 



Pa-Pb 
^P Av P Bp q a . 



pA* 



(2.15) 



°r = ~Qt- They form 



These vectors are orthogonal and satisfy Z 2 pa —Q 2 and X 2 pa Y 2 
a natural set of spacelike axes (with in the dilepton rest frame only spatial components). We 
will denote g M = q^ /Q, z» = Z^/y/^Z 2 , etcetera. Explicitly, one has in a collinear frame 



x 



y 



V2k 


V 


1 




V2 


• 

K 


1 


9. 



K q T 



K 



Vt 



0,0, 



Qi 



(2.16) 



where y l T = e^q^j. Note that since Z^ is a linear combination of the hadron momenta, 
it has in collinear frames no transverse components. The transverse vectors a^, thus, are 
orthogonal to Z. They are, in general, not orthogonal to q. One has, for example, q^ pa 
X M — (Qt/Q 2 ) <? m - Note that the second term is only order 1/Q suppressed. For an arbitrary 
four-vector a we define the perpendicular four-vector a_i_ as the projection of the transverse 
vector Or, using the projector 



ft 
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q^q" + z^z 



(2.17) 
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yielding 



a£ = gl a Tv = a T — g A 



Q 

Thus, any perpendicular vector satisfies 

o± • Q — a ± • z — 0. 



(2.18) 



(2.19) 



Note that X is the perpendicular projection of q. The vectors n + and n_, defined in 
Eq. ( |2.10D , can be expressed in terms of the set ( |2.16| ), 



1 Qt? 

n+ ~T2{ q + z ~^ x j 
n ~*72T~*~ ~q\ 



^2.20) 



Inserting these into the definitions of the projectors g± and g?, one derives the relation 

+ cfq T ) 



UV ^ , III/ 







from which one obtains for a general vector a: 



a i ~ air — 



Q 



2.21) 



(2.22) 



provided that ax ■ qr = 0(1). From this expression it is evident that for two arbitrary 
vectors a and 6, satisfying this condition, 



^2.23) 



Restricting oneself to leading order, the vectors and can be freely interchanged. In a 
higher order study, however, the difference will become important fTSfl . 

For the spin vectors the above definitions can be illustrated. In a collinear frame the 
spin vectors, satisfying Pa • Sa = Pb • Sb = 0, can be written as 



Sa 
S b 



A. 



M A 

2 P r 



A 



Ma 



'AT 



A 



P 



B 



B 



M, 



B 



A Mb S 



(2.24) 
(2.25) 



where A^, and A#, are the hadron helicities. The two-component vectors Sat, an d Sbt, 
give the transverse polarization. Since we consider pure spin states, they obey A 2 + = 1. 
For the spin vectors we have in a collinear frame St = [0,0, St]. The perpendicular spin 
vector is given by 



CM ^ CM 



St ' Qt 



Ol— t\,0[-),St + 0[ — 



kQ 



,Q 2 



(2.26) 
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where the longitudinal components follow from the transverse components by demanding 
Eq. fl2.19|) , and using Eq. Q2.16 ). If the spin vector would have been projected directly onto 
the XV-plane with g^ 1 , one would have got 



SfZS v 



[2.27) 



This differs from St in the transverse sector by 0(1), unless Qt = 0. 



C. Structure functions 



With the definition ( 2.15Q of a Cartesian set of vectors orthogonal to q, we can expand 
the lepton momenta in the following way: 

K = \(f + |Q(sin6»cos0 x M + sin # sin + cos 9 F), 

k% = |gf - |Q(sin 9 cos (px fl + sin 9 sin y^ + cos 9 2") . (2.28) 

Inserting these into Eq. (|2.2|) , and using some trivial goniometric relations and the com- 
pleteness relation g^ v = q^q v — z^z v — x^x u — y^y", we obtain 
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[I + cos 2 9) #f - 2 sin 2 9 z»z v + 2 sin 2 9 cos 20 {x^x u + \g^) 



+ sin^ 9 sin 20 x {ll y u} + sin 29 cos z iil x v] + sin 29 sin zS^y 



^2.29) 



where the symmetrization of indices, z^x v ^ = z^x" + z u x fM , is used. The six tensor com- 



binations in Eq. ( 2.29 ) are not only orthogonal to g, ensuring q^L^ = 0, but also to each 
other. 

It is convenient also to write the hadron tensor as a sum of products of tensors and scalar 
functions, called structure functions. From the properties of the electromagnetic current, 
one deduces for the hadronic tensor [Eq. (j2.3|)] the following conditions: 



q^W^ = 

W^(q;P A -S A ;P B -S B ) = W^(q;P A S A ;P B S B ) 
W^(q;P A S A ;P B S B ) = [W»"(q; P A S A ; P B S B )}* 



[Current conservation] 

[Hermiticity] 

[Parity] 

[Time reversal] 



(2.30) 



where = a M . The hermiticity condition, for instance, requires that the symmetric part of 
W^ v is real. In unpolarized scattering the constraints imply the expansion 



= -(Wo, - |W 2 , ) <?f + (W ,o + |W 2 ,o) z»z v - W 2A z^x^ - W 2 , 2 (x»x v + \gff) 

(2.31) 

where the four structure functions depend on the (four) independent scalars, or equivalently 
on Q, x A , x B , and Qt- Since we choose to work with the normalized vectors, the structure 
functions W 2A and W 2>2 contain kinematical zeros for X 2 « —Q T = of first and second 
order, respectively. In that they differ from the ones in RS || Eq. (2.5)]. To be precise: 



S 



our W21 is y/X 2 Z 2 times theirs, and our Wi,2 is —X 2 times theirs. The linear combinations 
multiplying —g^ and z^z v are often referred to as Wt and Wl, respectively. Inserting 
Eqs. ( |2.29| ) and (|2.31 ) into Eq. ( |2.1|) , one has for unpolarized Drell-Yan scattering 



da - 2 



2W 0)0 + W 2 , U - cos 2 9) + W / " 2 ,isin2^cos0 + iy 2 ,2|sin 2 ^cos20 . (2.32) 



d^qdVt 2s Q 

Due to the extra pseudo- vectors Sa and Sb, in polarized Drell-Yan there are several more 
structure functions. We will not give them in general. Later we will simply consider the 
ones that arise at leading order in 1 /Q in the cross section. 



III. FORMALISM 
A. The correlation function 

In order to calculate the hadron tensor, the diagrammatic expansion of EFP is used WA 



Each diagram is composed of soft nonlocal matrix elements, convoluted with a hard cut- 
amplitude. For example, in DY scattering the simplest diagrams (Born diagrams) are given 
in Fig. [| The nonlocal matrix element describing the nonperturbative long-range physics 
for a quark of flavor a inside hadron A, is given by the quark-quark correlation function 

(^/aUPaSa; k) = J e* k - x (PASA\^;\0)4 a) (x)\PAS A )c, (3.1) 

diagrammatically represented in Fig. ^. We will suppress the quark label a, the hadron 
label A, and the connectedness subscript c, whenever they are not explicitly needed. A 
contraction over color indices is implicit. More general, the correlation functions contain a 
number of quark and gluon fields, but at lowest order in 1/Q only the one above is relevant. 

A point that needs to be mentioned is the color gauge invariance of the correlation 
function defined above. For a proper gauge-invariant definition of a nonlocal matrix element 
a color link operator, L(0,x) = V exp[ig f$ ds ■ A(s)], must be inserted. As will be seen 
below, we only need the integral of $ over k~ , i.e., the nonlocality is restricted to the plane 
x + = 0. Going further and integrating also over transverse momenta, one is only sensitive 
to the nonlocality in a lightlike direction (say x~). With an appropriate choice of path (a 
straight link) and lightcone gauge (A + = 0), the link operator becomes just unity. When 
transverse momentum is observed, this no longer is the case, since one becomes sensitive 
to transverse separation X?, although still x + = 0. Fixing the residual gauge freedom that 
affects At can be achieved by imposing boundary conditions [[0J. A path from to x can 
then be constructed such that the link operator becomes unity after gauge-fixing. How the 
choices of gauge and path connect to the proof of factorization remains an open question. 



B. The Dirac structure of the correlation function 

In order to analyze the diagrams in Fig. ^| for DY scattering, we need to investigate the 
Dirac structure of the correlation function. This can be done by making an expansion in 
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an appropriate basis. Constraints on the correlation function come from hermiticity, parity 
invariance, and time reversal invariance, 

&(PS; k) = 7° $(PS;J;)_7 [Hermiticity] 

$(PS; k) = 7° $(P -Sj k) 7° [Parity] (3.2) 

$*(PS; k) = 75 C $(PS; k) &j 5 [Time reversal] 

where the charge conjugation matrix C = 27 2 7°, and fc M = k^. Choosing the Dirac matrix 
basis 1, ijs, 7 M , 7^75, and ia^^s (note that = 7°r7°), the most general structure 
satisfying these constraints is 

$(PS; k) = At 1 + A 2 P + A 3 j& 

+A Al5 p + A 5l5 [P, t$] + A 6l5 [k', ft 

+A 7 k • 5 75 P + A s k ■ S lb U* + A 9 k ■ S l5 [P, U]. (3.3) 

Hermiticity requires all the amplitudes Ai = Ai(k ■ P, k 2 ) to be real. Note the presence of 
the amplitude A 9 which is left out in Eq. (3.4) of Ref. 0. 

The basic assumption made for the correlation function is that in the hadron rest frame 
the quark momentum k is restricted to a hadronic scale A, explicitly k 2 and k ■ P are of 
0(A 2 ). In a frame where the hadron has no transverse momentum, the momentum k is 
written as 



k ~~\~ krjp 

2xP+ ' 



xP + ,k T 



(3.4) 



with the lightcone momentum fraction x = k + /P + . The restrictions on k 2 and k ■ P imply 
that also k 2 -, = —k 2 + 2xk ■ P — x 2 M 2 is of 0(A 2 ). Considering diagram 0a, one sees 
easily that momentum conservation on the hard vertex implies q~ = k~ + k^. However, 
k^ ~ Pg, whereas k~ ~ M\jP\^ which is down by a factor ~ M 2 /Q 2 in any collinear frame. 
Therefore, for hadron A, one is led to study / dk~ $(P5*; k), or equivalently its projections 
/ d/c _ Tr[r$]. These latter quantities do not carry Dirac indices anymore, but because of 
the T-matrices, they do have a specific Lorentz tensor character. Defining the projections 

$[r](z,fc T ) = - J dk~ Tr [r $] (3.5) 

= \l ^T^i exp[i(xP+a;- - k T ■ x T )\ (P5|^(0)IV(0, x~, x T )\PS), (3.6) 

l J ZTT (27TJ 

one has for instance the vector projection 

$[ 7 +] = J d(2k ■ P)dk 2 5 (k 2 T + k 2 - 2xk ■ P + x 2 M 2 ) (A 2 + xA 3 ), (3.7) 

which is 0(1). Other projections, e.g., the scalar 

$[1] = _ J d{2k ■ P)dk 2 6 (k 2 T + k 2 - 2xk ■ P + x 2 M 2 ) A 1} (3.8) 

contain an integral of 0(1) multiplied by a factor 1/P + . In the cross section, this factor will 
give rise to a suppression of order 1/Q. In this way it is seen that the leading contributions 
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come from the Dirac structure where the number of +-components minus the number of 
— components is largest (that is, 1). They are parametrized as 

<% + 7 5 ] = 9il(x, k 2 T )X + g 1T (x, k 2 T )^j^, (3.9) 



$[ia i+ l5 ) = h 1T (x,k T )S T + 



hi L (x, k T )\ + h 1T (x, k% 



kx ■ St 

M 



defining six real distribution functions per flavor, depending on x and k\. These encode the 
leading behavior of the quark correlation function. An expansion can be made in terms of 
local operators which have twist two and higher. 

In the diagrammatic expansion for the DY hadron tensor (with A + = for the lower 
blob), correlations will appear with, in addition to the quark fields, transverse gluon fields 
in the soft part, ip(0)A T (y)ip(x), etc.. They can be analyzed in the same way, and turn out 
to be suppressed by one order of 1/Q for each additional gluon field. 

In summary, for a leading-order DY calculation, one is to use 

J dk~<5> = |$[ 7 +] 7 - + |$[ 7 + 75] 757" + !$K + 7 5 ] z 75( 7l + . . . , (3.10) 

where i is a transverse index (i.e., i = 1,2), and the dots represent projections that will 
come in only at 0(1/ Q). At the leading-order level, one is left with only three projections, 
which have vector, axial-vector, and axial-tensor character, respectively, and which can be 
parametrized by six distribution functions. 



C. Antiquarks 



The antiquark correlation function, describing the antiquarks of flavor a, is given by 
(contracting over color indices) 



d A x 



ik-x 



(p A s A \4 a) mf\x)\p A s A ), 



(3.11) 



Also the antiquark momentum k can be written as in Eq. ( |3.4[) . Its Dirac structure can be 
analyzed likewise. We define the antiquark projections 



$\T](x, k T ) 



dk Tr 



r $ 



(3.12) 



1 r dx d 2 XT 



.., exp[i(xP + x- - k T ■ x T )} (PS\Tr ^(0)^(0, x~,Xr) \PS). (3.13) 
Using the charge conjugation properties of Dirac fields and hadron states, we deduce 

^- a/A = -C- l (^ a/A ) T C. (3.14) 

Upon demanding charge conjugation invariance of the distribution functions, i.e., the quark 
distributions in the antihadron A are the same as the corresponding antiquark distributions 
in A, we obtain the expressions 
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$[ 7 +] = A0r,fc!), 

$[7 + 75] = ~9il(x, k^)X - g 1T (x, k 2 T )- 7 



<S>[ia i+ j 5 ] = h 1T {x,k 2 T )S T + 



M ' 



kj< ■ St 



M 



iXjrp 

~M' 



(3.15) 
(3.16) 

(3.17) 



We note that the anticommutation relations for fermions can be used to obtain the 
symmetry relation 



$ij(PS; k) = -^(PS; -k). 
For the distribution functions this gives the symmetry relation 

fi(x, k 2 T ) = -fi(-x, k 2 T ), 
and identically for g 1T , h XT and h± T , whereas 

§il(x, k 2 T ) = g 1L (- x , k T ), 



(3.18) 



(3.19) 



(3.20) 



and identically for h\ L 

Finally, we note that hadron B can be treated in the same fashion. However, since we 
have chosen to work in the collinear frames where the third axis lies opposite to the direction 
of hadron B, the role of the +- and — components for B must be interchanged as compared 
to A. So, in Fig. |2]a the upper blob reduces to / dk£§(PBSB', kh), for which one has to use 
the gauge A~ = 0. 



IV. FREE-QUARK ENSEMBLE 



The leading fc T -integrated distributions fi(x), g%(x), and hi(x), have a parton model 
interpretation as the longitudinal momentum, helicity, and transversity distribution, respec- 
tively. In this section we show how, for a free-quark ensemble, this identification can be 
generalized to the fc^-dependent distributions. 

It is instructive to calculate the correlation function for a free-quark target of flavor a. 
This is given by 



(®a/a)ij(p s; k) = 5 4 (k - p) Ui(p, s)uj{p, s) = 5 4 (k - p) 
where the momentum and spin vector are parametrized as 
k 



+ m) 



(4.1) 



ij 



m2 + k ? k\k T 



2k+ 

s = A a n fc + s at = X a 



k T — m 2 k + k T 

2m k + ' m ' m 



+ 



k T ■ s aT 

k+ 



, 0, SaT 



(4.2) 
(4.3) 



We identify the lightcone helicity A a , helicity vector n^, and transverse polarization s at 
(transverse in the sense that s at ■ k — s at ■ — 0). Note that the lightcone helicity vector 
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rik, satisfying rik ■ k = and n\ = — 1, acquires its conventional meaning Eq. (2-49)], 
either if = Ot, or in the infinite-momentum-limit k + — > oo with fc^ fixed. One checks 
that X 2 , + s 2 T = —s 2 = 1. With the simple form (fOJ) it is straightforward to calculate the 
leading projections for a free-quark target with non-zero transverse momentum, 

l - J dk~ Tr[ 7 + $(p s; k)] = 5 - 6 2 (k T - p T ), 

i | dAT Tr[ 7 + 7 5$(P s; k)] = 5 ( ^ - lj <5 2 (fc T - p r ) A a , (4.4) 



i | dfc" Tr[i<7 i+ 7 5$(p «; fc)] = M ^ ~ 1 J * 2 (fc r - Pt) s{ t . 



This simple example of a quark target can be generalized to the case in which the 
hadron is considered as a beam of non-interacting partons of total momentum P and angular 
momentum S. This is tantamount to inserting free-field plane-wave expansions into the 
correlation function ( |3.1| ). One gets (summing over a, [3 — 1,2) 

<S>ij(PS; k) = 25{k 2 - m 2 ) [9{k + )uf\k)Vp a {k)uf\k) - e{-k + )vf\-k)Vp a {-k)vf\-k) . 

(4.5) 



The functions V and V are given by 

1 r dr' d 2 V 

Vp a {k) = V 0a (x, k T ) = _ J ^-^X (PS\bl(k'Mk)\PS), (4.6) 



The form of quantization one can choose to be instant-front, as well as light-front quantiza- 
tion, since for free fields they are equivalent JO! . As for the choice of coordinates, the use 



of lightcone coordinates is convenient, because of the integration over k that is needed in 
deep inelastic processes. The Dirac structure can be parametrized as 

u^\k)V Pa {k)u^\k) = V{k)U& + m) ( 1 + 7 2 5/(fc) ) , (4.8) 
vW(k)Vpa{k)it«\k) = V{m - m) (^ 1 +7 2 57(fc) ) , (4.9) 

in terms of positive definite quark and antiquark probability densities V(k) and V(k), and 
spin vectors s M (/c) and s M (/c) . Inserting the free-field expansion in the current expecta- 
tion value (PS'|?/'(0)7^(0)|PS') = 2P fl (N - N), where N and N are the total number of 
quarks and antiquarks, respectively, one obtains from the +-component the normalizations 
Jq 1 dx J d 2 kTV(x, fc|) = N and Jq dx J d 2 kTV(x, k\) = N. The average quark spin vector 
s M (fc) is parametrized by the helicity density X a {x, kr) and transverse polarization density 
s a T{x, kx) by expanding s M as in Eq. (|4.3|) . A similar parametrization is used for ~s^{k) in 
terms of A a (x, k T ) and s sT (x, k T ). 
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Integrating Eq. (|4.5|) over k one obtains the result for a free-quark ensemble, 



2k+ w ' \ 2 
-x, 



- e(-x) V (-*f T \ l6 + m) ( . (4.10) 



This gives (for x > 0) 

dk-Tr[ 7 + ^PS;k)]=P(x,k 2 T ), 

X - J dk~ Tr[ 7 + 75 $(P£; k)] = V(x, k 2 T ) X a (x, k T ), (4.11) 

l - J dk- Tr[i<r i+ l5 ${PS; k)} = V(x, k 2 T ) 4 T (x, k T ). 

A comparison with Eq. ( |3.9| ) yields 
V(x,k 2 T ) = f(x,k 2 T ), 

V(x, k 2 ) X a (x, k T ) = g 1L {x, k 2 T ) A + g 1T (x, k 2 T ) kr 'J T , (4.12) 



V(x, k T ) si T (x, k T ) = h 1T (x, k 2 T ) S T + 



.2,1 / . 1 2 \ ■ St 



M 



~M 



which shows how for x > the functions giL, #it, ^itj h± L , and /i^ r , are to be inter- 
preted as quark longitudinal and transverse polarization distributions. The function hj; T 
was omitted in the paper of RS. Specifically for non-zero transverse momenta of the quarks, 
it becomes relevant. For the antiquarks the same relations hold between the antiquark 
probability density V, helicity density A„, and transverse polarization density s^t, on the 
one hand, and the antiquark distributions on the other hand. Extending to all x, re- 
sults are obtained in accordance with the symmetry relations in the previous section, e.g., 
f(x,k T ) = 9(x)V(x,k T ) — 6{—x)V{—x,k' T ). For on-shell quarks, relations exist between 
the distributions due to constraints from Lorentz invariance ||14|| . The scalar V can only 



depend on 2k ■ P = m 2 /x + M 2 x + k^/x and the pseudovector can be parametrized by 
two functions depending on the same combination of x and k\. 

We note that for the leading-order matrix elements the free-field results can be used 
to provide a parton interpretation, even in the interacting theory, because the distribution 
functions can be expressed as densities for specific projections of the so-called 'good' com- 
ponents of the quark field; ip + = A + ip, where A + = |7~7 + . In light-front quantization a 
Fourier expansion for the good components (at x + = 0) can be written down in which the 
Fourier coefficients can be interpreted as particle and antiparticle creation and annihilation 



operators The different polarization distributions involve projection operators that 

commute with A + 

At subleading order, the analysis of the quark-quark correlation functions leads to a 
number of new distribution functions. For free quarks, they can also be expressed in the 
quark densities and thus they can be related to the leading distribution functions. However, 
it turns out that the presence of nonvanishing quark- quark-gluon correlation functions causes 



deviations from the free-field results [18 
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V. RESULTS 



A. Hadron tensor 



Using the EFP-expansion, the leading-order Drell-Yan hadron tensor in the deep inelastic 
limit, and with Qt = O(A), can be written as the sum of the quark and antiquark Born 
diagrams in Fig. |2|. First, we will calculate diagram [^a, in which a quark of hadron A 
annihilates an antiquark of B. It reads 

WTuark = \ E h- a e 2 a J d A k a d A k b 5\k a + k b -q) Tr [$ a/A (P A S A ; k a ) 7" $ b/B (P B S B ; k b ) 7 

a, b 

(5.1) 

where a (6) runs over all quark (antiquark) flavors, and e a is the quark charge in units of e. 
The factor 1/3 comes from the fact that the quark fields in both the correlation functions are 
traced over a color identity operator, which is appropriate since only color-singlet operators 
can give non-zero matrix elements between (color-singlet) hadron states (Wigner-Eckart 
theorem). However, since the diagrams we consider have only one quark loop, one has 
only one color summation, leading to a color factor 1/3. Using the boundedness of quark 
momenta in hadrons as discussed before, one finds k+ ^> k£ and k^ ^> k~ . Thus, the delta 
function can be approximated by 

5\k a + k b -q)n 6(k+ - q + )5(k b - - q-)5\k aT + k bT - q T ). (5.2) 

It follows that indeed we are sensitive to J dk~^ A (k) and / dk + & B (k). Furthermore, the 
trace in Eq. ( |5.1|) can be factorized by means of the Fierz decomposition 



ljihk + (iy 5 )ji(ij5)ik - (l a )ji(la)ik - (l a/ y5)ji(lalb)ik + \{io ^75) ji(ia^ 
Hi { %i{i v} )i* + (tNO^YSOi* + {i<J a{ ^)^° u} ^)ik + • • • , (5.3) 



where the dots denote structures antisymmetric under the exchange of \i and v. We keep 
only the leading projections, as they where found in Sec. [Til]. This leads to 

^quark = ~\ E J d 2 k aT d 2 k bT 5 2 (k aT + k bT - q T ) 

a,b 

§a/Ab + }§b/B[T] + ^ a /Ah + l^b/Bb^l^j9T 

+$a M k i+ 7 5 ] ® b/B [ia j - l5 ] (jnVi - m&) + O , (5.4) 

where the ^/^-projections are taken at x = x A and k a T, the ^/^-projections at x = x B and 
k bT . The next step would be to insert Eq. ( |3.9|) and their hadron- B antiquark counterparts. 
In this leading-order calculation we may interchange the subscript T by a _L where we wish, 
because of Eqs. ( |2.21| ), ( |2.22|) , and (|2.23| ). The resulting expression contains convolutions 
like 
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J d 2 k aT d 2 k bT 5 2 (k aT + k bT - q T )F(k 2 aTl k 2 bT )k% ± . 



(5.5) 



They are not of the desired form, since the (perpendicular) Lorentz index is carried by a 
convolution variable. In order to write the hadron tensor in terms of structure functions, the 
Lorentz indices must be carried by external vectors, like X M . In the appendix we describe 
the method we used to project perpendicular Lorentz tensors, like (|5.5| ), onto the XF-basis. 



For instance, (5J5) can be written as (up to order 1/Q 



d 2 k aT d 2 k bT 5 2 {k aT + k bT - q T )F(k 2 aT , k 2 bT ) 



bT 



2Qi 



(5.6) 



Before presenting the results, we mention that the antiquark diagram ^ can be obtained 
from the quark diagram |2|a by making the replacements^ a «-> a, f «-> /, g «-> — g, h <-> h, 
where g is generic for the quark axial- vector distributions, etcetera. It is convenient to define 
the convolution of two arbitrary distributions dx(x, fey) and d~2(x,k T ) (possibly multiplied 
by an overall function of k 2 aT and k 2 T ) 



I[d x d 2 ] = ~^8 m e 2 a J d 2 k aT d 2 k bT 5 2 (k aT + k bT - q T )di(x A , k 2 aT )d 2 (x B , k 2 bT ) 

a,b 



(5.7) 



summing over quark and antiquark flavors, with the prescription that if its argument is of 
antiquark (quark) nature, then d\ (d 2 ) is to be replaced, either by d\ (d 2 ) if it concerns a 
vector or axial-tensor distribution, or by —d\ (—d 2 ) if it concerns an axial- vector distribution. 
We split up the general leading-order hadron tensor W J,l '(S A , S b ) in four cases, corresponding 
to unpolarized, longitudinal-longitudinal, transverse-longitudinal, and transverse-transverse 
scattering, 



^(0,0) 
W^(X A ,X B ) 

W^{S AT ,X B ) 



-W T g?, (5.8) 

- {Wt + \V T LL X A X B ) of - \V^X A X B {x^ + ±af ), (5.9) 

- W T - V? L ( X ■ S A± )X B ] of + V 2 T 2 L (x ■ S AA _)X B (x»x v + |of ) 



—U 2 2 X B 



W» V (S, 



AT, &BT 



[x ■ SAA.)g±\ , 

r2 TT i 



(5.10) 



W T - Vl TT (S A± ■ S BX ) + V^ T (x ■ S AX )(x ■ S BX ) 



9? 



+ 



^2L2 1 (S A ± ■ S B ±j - . 2 2 



ITT 

-ATT 
2,2 

■BTT 

'2,2 

T TT 



V 2 r{x-S AA _){x-S B± ) (x^ + hT) 



+U^ 1 (x-S B± ) 



x^S u A } ± - (x ■ S AX )g^ 



+l T ?o'Ux.S A± )[x^S$ ± -(x-S B± )g d a '] 
-U 22 [Saj-Stsx ~ (Sa± ■ S B ±)g^ , 
where the structure functions, depending on x A , x B , and Qt, are given by 



(5.11) 



3 It is seen from Eq. ( |5.4| ) that this amounts to replacing A <-> B, rendering a result, symmetric 
under the exchange of the two hadrons. 
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W T = I[hhi (5.12a) 
V^ L = -AI[g 1L g 1L ], (5.12b) 



Q 2 T ) M A M B 



V 2 L 2 L = I[ a + P- (5.12c) 



Ql J M\MbQt 



2M a Q t 

VL L = /[( PQ T + a> + aP- 2/3 2 - ^LZ®1) -^-]. (0.12c) 



2M B Q T 

+ ma -e)-2(a>-P ) + ^) I *ff^ ] , (5-120 



^-nl-Qi^^-^f)^, (5, 2g) 



V2ir_j„ 2 , , 2 2(a + /3)(q-/3) 2 ( Q -/3) 4 \ fej- r fej- r 



q 2 ^ / m\m% 



t/1 TT i j j ATT | jtBTT 
V 2,2 ~r ^2,2 ~r ^2,2 



u(n*( n +R\ ,,2 , 5(a + p)(g-py 2(q-/3) 4 \ hj T hj T 
+ {Qt(* + P)-^ -4/3 + ^ — I——], (5 ' 12j) 



tATT ttBTT T J n 2 o„ | ^IT^jr / n 2 i («-fi) 2 ^ feujlT 



+ [«( - /3) - 2a * + ^ + ^)^] l (6.12k) 



rrTT , ItjATT , ItjBTT _ j\ h T h 1T h 1T Jlyrh 1T 

^2,2+2^2,2 +2^2,2 ~ J [^lT^lT + « ^2 + /* 2 M| 

_(,„_„?_ 2(a + /3)(a-/3) 2 (a-/3) 4 \ 

^ ^ j Q£ J8M 2 M 2J ' 

where a = fc 2 T , and /3 = fe 2 T . In the last three equations we only gave the combination 
of structure functions that will show up independently in the cross section, to be discussed 
in the next subsection. As for the nomenclature, we based ourselves on the RS-notation, 
except that we always denote both hadron polarizations by the TL-superindices. So, for 
instance, their V£ corresponds to our V^jf , and their Z7 2 ,2 equals our U^J- The subindex T 
stands for the linear combination (0, 0) — |(2, 0). Note that the hadron tensors (|5.8|) - (|5.11 ) 
manifestly conserve the electromagnetic current, i.e., q^W^ = 0. 
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B. Cross sections 



It is now a straightforward matter to contract the hadron tensor with the angle-dependent 
leptonic tensor ( [2.29 ). For the four combinations of polarization we deduce the following 
leading-order differential cross sections: 



da(0,0) 


a 2 


d^qdVL 


2sQ 2 


da(X A , Xb) 


a 2 


d^qdtt 


2sQ 2 


da(S AT , Ab) 


a 2 


d 4 qdn 


2sQ 2 


da(S A r, Sbt) 


a 2 


d^qdVl 


2sQ 2 



W T (1 + cos 2 



w T + \v^ L \ A \ B 



[l + cos 2 6) + \V^\ A \ B sin 2 6 cos ! 



(W T + V^ l A b cos0a) (1 + cos 2 6») + %v£f\ B coa<i>ABm 2 8co8 
+U 22 \ B sin 2 9 cos(20 - A 

(W T + Vt TT cos((J) A - 4>b) + V^ tt cos0 a cos0b)(1 + cos 

+ ^V 22 cos(p A cos 0# sin 9 cos 2(f) 
+\{ V 2,2 T + U 2,2 TT + U 2,2 TT ) cos (<^ - 4>b) sin 2 9 COS 2(f) 
+U U 2,2 TT ~ U 2,I T ) sin (0A - 4>b) sin 2 9 sin 20 



(5.13) 
(5.14) 

2(f) 

(5.15) 



HU'2,2 + M? T + \ U 2,2 rT ) sin " 9 cos ( 2 <^ " ^ - 0b) 



1 rrB TT 



(5.16) 



where A (0s) is the azimuthal angle of S^x (Sb±), i-e., cos0 A = — x ■ S A ±_ ~ • Sat/Qt 
for pure transverse polarization (|iSUt| = 1)- Eq. (|5.16|) shows that the structure functions 
V 2 1 2 TT , U 22 TT , U 22 TT , an d cannot separately be extracted from the experiment, but 
only in three particular combinations. This is due to the fact that the corresponding tensor 
structures in Eq. ( 5.11 ) are not independent, because of the relation 

g p °(2x^ + fff) - x"(x^gt - x°g^) - x°(x^g? p - Pg?) - (g^g^ ~ fig?) = 0, 

(5.17) 

which can be proven by using g^ = —x^x u — y fJ "y v '. One easily checks that the angular 
functions of Eqs. (|5.13| )- (|5TT6|) have no further dependencies. 

In order to circumvent normalization problems, in spin experiments one usually considers 
the asymmetries 



A? 



a(S A , S B ) - cr(S A , -S B ) 



SaSb ct(S a ,Sb)+ct(S a ,-Sb)' 
For example, the longitudinal- longitudinal asymmetry follows from Eq. (|5.14 ), 



(5.18) 



.4 



-^A^B 



V^ L sin 2 9 cos 20 



W 7 



1 + cos 2 9 W 7 



(5.19) 



where the structure functions can be found in Eqs. ( |5.12a| ), ( |5.12b ), and ( 5.12c| ). The 
angle-independent term is a generalization of the qr^-integrated longitudinal-longitudinal 
asymmetry [^T] . 
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VI. DISCUSSION AND CONCLUSION 



In this last section we will look more closely into the results we derived, first by comparing 
the cross sections with those integrated over q T , secondly by considering the limit Qt — > 0, 
and finally by considering a Gaussian fc^-dependence of the distribution functions. 



A. Qy-integrated results 

If the Qr-dependence is eliminated by integrating over the transverse momentum of 
the produced lepton pair, one recovers the lightcone momentum, helicity, and transversity 
distributions of Eqs. (|T|), flip , and flLg); 

fi(x) = J 'd 2 k T h{x,k 2 T ), (6.1) 

9i( x ) = / d2k T 9il(x, fey), (6.2) 

(6.3) 



hi(x) = J d 2 kq 



k 2 

h 1T (x, k 2 T ) + T^^rrO, k 2 T ) 



and similar antiquark distributions. This can be seen in two ways. It is easiest to return 
to Eq. ( |5.4|) and observe that the delta function is absorbed by the qr-integration, so that 
the transverse intergrations over k a x and fc^r separate. One is left with three independent 
structure functions, 

W T (x A , x B ) = § e l h{xA)hix B ), (6.4a) 

a 

V T L (x A ,x B ) = ~^Yl e l 9i(xa)9i(xb), (6.4b) 

a 

U^(x A , xb) = \Y1 e a h i{x A )hi{x B ), (6.4c) 



with the by now well-known prescription that if a is an antiflavor, the vector and axial- 
tensor distributions have to be replaced by their charge conjugated partners, the axial- vector 
distributions by minus their anti-partners. The names of the structure functions in Eq. ( |6.4|) 
refer to the Lorentz tensor structure they multiply in the hadron tensor (and hence the 
angular distribution in the cross section). That is, Wt in the hadron tensor is the coefficient 
of —g 1 ^, etcetera. The same expressions are obtained if one integrates the full results in 
Eqs. (|5.8|) - (|5.11|) over qT- It is then seen that Wt and V T are the integrals of Wt and V^ L , 



respectively, but that U % \ originates from the combination + ^U^ T + |^2 TT + \V%2 T - 
Contraction of the hadron tensors with the lepton tensor leads to the following leading- 
order four-fold differential cross sections ||: 

J^ML = ^w T (i + ^e), (6.5) 



dx A dx B dQ AQ 
da(\ A ,X B ) a 2 — itjLL 



dx A dx B dQ 4Q 5 



(W T + \V T Xa\b) (1 + cos 2 6) , (6.6) 
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da(S A Ti Ag) 
dxAdxsdQ 
dcr(S A T, Sbt) 
dxAdxsdVt 



a 

W 

a 2 

4Q 2 



W T (1 + cos 2 



M + cos 2 



+ u 2 2 sin 



6* cos(20 — 4>a — 4>i 



(6.7) 
(6.8) 



Note that no absolute azimuthal angles occur, but only the relative angles <fi — 4>a and 
— 4>b- Two remarks are in place here. Comparing with the full cross sections ( |5.13| - |57L6"1) , 
we observe the disappearance of the sin 2 6 cos 20-term in the longitudinal-longitudinal cross 
section, and hence in the asymmetry ( |5.19 ). Secondly, we see no polarization-dependence in 
the integrated transverse-longitudinal cross section, and hence no asymmetry. Indeed this 
asymmetry is suppressed by a factor ~ 1/Q M. 



B. Q 3 



Only a few of the structure functions are nonvanishing at Qt = 0, the other ones having 
kinematical zeros (see also [JL5|J). Therefore, it is worthwhile to consider the limit Qt — > of 
the structure functions in Eq. ( |5.12| ). In general we can do this, rewriting the convolutions 
in Eq. (|5.7| ) by making a transformation to the momenta k T = \{k a r — kbr) and Kt = 

k a T + k bT , 

I[dJ 2 ] = |E e ' / d 2 k T dJxA, (k T + \q T f) d 2 (x B , (k T - \q T f) (6.9) 

= IE e « / tfkrfaixA,^) d 2 (xB,J4) + O(Q^) = I [d 1 d 2 ] + O{Q^), (6.10) 

assuming the distribution functions to be sufficiently well-behaved to justify the Taylor 
expansion of the integrand. For instance, for V 2 L 2 L of Eq. ( |5.12c| ) one finds in the limit 
Qt —* (omitting the flavor sum) 

' (q T -k T ) 2 



V 2 L 2 L <X 



d 2 kn 



qWt 



2 k T 



Qt 



+ \Ql hi L (x A , (k T + \q T f) ht L (x B , (k T - \q T )' 



Qt 



d 2 k T - k Ti k Tj ) h± L (x A , k 2 T ) hi L (x B , k 2 T ) + 0(Q 2 T ) = 0{Q 2 T ). (6.11) 



So V 2 2 has a kinematical zero of second order, which is the natural behavior, because it 
multiplies a tensor that is quadratic in x and y. In the same fashion all structure functions 
can be treated. Since this is a rather cumbersome procedure, we will illustrate below the 
behavior by considering a Gaussian k\- dependence, from which the order of the kinematical 
zeros is simply read off. We obtain only four structure functions without a kinematical zero, 
in agreement with the results of RS, 



t\q t =o 



Vr, 



LL 



Qt=0 



Vr, 



ITT 



Wi/i], 

= -4Io\9il9il], 
- MA- 2 9lT ^ 1T 



2.1/ ,.!//; ' 



■JJ TT _l UjATT , ltj 
u 2,2 ' 2 2,2 ' i u ' 



IttB TT\ 
2 U 2,2 J 



?T=0 



Io[{hiT + ^hi T )(h lT + k 



2M|^ t) 



(6.12a) 
(6.12b) 

(6.12c) 
(6.12d) 
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An easier way to obtain this result is to start from Eq. ( |5.4j) , and put Qt = from there. 
One then never picks up the other structure functions in the first place. The cross sections at 
Qt = can be obtained by insertion of these structure functions into the explicit expressions 
given for the cross sections in Eqs. (|5.13|) - (|5.16|) . 



C. Gaussian transverse momentum distributions 



It is instructive to consider a Gaussian fc r -dependence, 

d{x, k T ) = d(x, 0) exp (-r 2 k T ) , (6.13) 



where the transverse radius r in principle depends on both the particular distribution func- 
tion, and on x. One can explicitly perform the convolution integral ( |5.7|) (for simplicity we 
will omit the color factor and flavor summation), 

/[dida] = - 71 exp \-Q 2 T ^ B 2 ) d x {x A , 0)d 2 (x B , 0), (6.14) 
n + r B \ r A + r Bj 

being regular in Qt = 0. One can identify an 'average' transverse size r, given by r -2 = 
r A 2 + r B 2 - We find for the structure functions, 

W T = I[hh], (6.15a) 
V^ L = -AI[g 1L g 1L ], (6.15b) 



v »=Ss^S? 11 ^ (6 ' 15c) 



^ L = -w A ht^^I- (6 - 15d) 



2^,2 



r2 TT _ 2Q T r A r B 

" M\M\ (r 2 + r 2 ) 



Vll 1 = ^7T7^rZ^nhf T hU (6.15i) 



t/1 TT , j j ATT i rrBTT _ r A ttt 7i i , ' i/ /T/.-L /, 

',2 + ^ 2)2 " ^ (r 2 + r 2 )2 + J^ (r 2 A + r 2 )2 7 M^' 



* ^-^^-^A], (6.15J) 



2M 2 M 2 (r 2 + r|)3 

fMTT rrBTT _ r A rr. 7_L i Qt ' Li in I /, 
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Q 



2M\M 2 B 



[ri + r 



ttTT i IttATT , ItjBTT 
U 2,2 2 2,2 2 2,2 



2M 2 ( 



1 + 
1 



2M 2 A M 2 B 



{r\ + r 




(6.15k) 



(6.151) 



These explicit results illustrate the kinematical zeros, and can be used to obtain their order. 
For instance, Eq. ( |6.15c| ) is an illustration of the result in Eq. (|6.11|) . The expressions can 



also be used to illustrate the expected behavior of the longitudinal-longitudinal asymmetry 
in Eq. ( |5.19| ), which becomes 



A 



I[giLgiL\ . sin 2 6 cos 20 Q\ 



+ 



r\r 2 B 



I[h{ L h 



ill 



I[hh} ' 1 + cos 2 6 M^M B {r\ + r 2 B ) 2 I[hf x \ ' ( ° 6) 

The second (angle- dependent) term in the asymmetry starts off with Q T , and is proportional 
to the average transverse radius squared of the functions hj~ L and h\ L . As a second example 
we note the structure function V^ L , being directly proportional to the transverse radius 
squared of the helicity distribution g 1L for hadron B. 

Finally, we mention the CERN NA10 experiment [22|, which indicates a possible 
sin 2 6 cos 20-asymmetry in unpolarized DY scattering at measured Qt, not suppressed by 
powers of 1/Q. Such an asymmetry, however, does not appear in the leading-order re- 
sult ( |5.13| ). It is, however, important to note that any kind of polarization in the beams 
leads to such an asymmetry. Recently, Brandenburg, Nachtmann, and Mirkes [23] have sug- 
gested that the nontrivial QCD vacuum structure could be responsible for the asymmetry. 
Another suggestion is that higher-twist effects are responsible p4| . 



D. Conclusion 

In this paper we have investigated the leading-order polarized Drell-Yan cross section 
at measured transverse momentum of the lepton pair Qt A. It is sensitive to the in- 
trinsic transverse momentum of the quarks and antiquarks in the colliding hadrons. The 
results can be written in terms of six quark and six antiquark distributions, depending on 
longitudinal lightcone momentum fraction x, and transverse momentum k T . For each flavor 
they represent the lightcone momentum, helicity, and transverse polarization distributions. 
Measurements of these distributions require a study of the asymmetries in doubly-polarized 
DY scattering, using longitudinally and transversely polarized beams. 

The quark distributions of Sec. [I IT] incorporate the intrinsic fc^-dependence. Although 
this involves higher-twist operators, the presence of a second scale Qt in the process en- 
ables one to find observable effects that are not suppressed by the large scale Q. Similarly, 
fer-dependent correlation functions can be used to analyze other deep inelastic processes, 
particularly semi- inclusive DIS. They will then appear in convolutions with the quark frag- 
mentation functions. At this stage factorization has not been proven, so we do not know if 
the distribution functions are universal. 
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APPENDIX A: PERPENDICULAR PROJECTION METHOD 

Consider the following convolutions: 

J[fcM • • • Kl\ = J d 2 k aT d 2 k bT 5 2 (k aT + k bT - q T )F(k 2 aT , k 2 bT )k^l . . . kfc. (Al) 

The perpendicular vectors ki± are taken from the set {k a ±, k b ±}. In practise we only needed 
the cases where n = 0,1,2,3, and n = 4 with two of the indices symmetrized. We want 
to project these perpendicular tensors onto a basis of Xy-tensors, multiplied with scalar 
functions of Qt- 

Take the simplest non-trivial case of n — 1, 

AKA = J d 2 k aT d 2 k bT 5 2 (k aT + k bT - q T )F(k 2 aT , k 2 bT )k^. (A2) 

Since the only perpendicular direction available is ~ X M , this expression must be pro- 
portional to x^. The proportionality constant is easily obtained by contracting with x^, and 
using x 2 = —1. We get (discarding corrections of 0{1/Q 2 ) in this appendix) 

AKA = -^Jfi ■ k iA = ^AQt ■ k 1T ]/Q T , (A3) 

making use of Eq. (|2.23|) . We choose to write the integrand as a function of Qt, k 2 aT) and 
k 2 T , by making use of the relations 

q T -k aT = \ (Q 2 T + k 2 aT - k 2 bT ) , ( A4) 

Qt 'kvr=\ (Q 2 T - k 2 aT + k 2 bT ) , ( A5) 

k aT -k bT = \ {Q 2 T - k 2 aT -k 2 bT ). ( A6) 

For the more intricate cases, n > 1, the procedure is essentially the same. Instead of £c M , 
and y >1 , we will use the more convenient building blocks £ M , and g*A = —x^x u — y^y u . For 
n = 2 we find 

JK^A = ( z ) ( \\ ) ( ' kM ) (A7) 

= (2x»x» + gT)A(QT ■ k 1T )(q T ■ k 2T )]/Q 2 T - {x»x v + g^) J[k 1T ■ k 2T ], (A8) 

where the matrix notation of the first line should be clear from the second. In the same 
notation, the n = 3-case reads 
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J[k 1± k 



v k p 1 



rfi^rf^rfP ~y* 

iAj iAj tAJ y %KJ 



x u g^, x p g u £ 



/4 1 1 1\ / J[(q T ■ k 1T )(q T ■ k 2T )(q T ■ k 3T )}/Q 3 T \ 



110 
10 10 
V 1 1 J 



J[(q T ■ k 3T )(k 1T ■ k 2T )]/QT 
-J[(q T ■ k 2T )(kiT ■ k 3T )]/Q T 
V -J[{Qt ■ k 1T )(k 2T ■ k 3T )]/Q T 



(A9) 



Finally, for n = 4 and two indices symmetrized, we find 



i i _ - \ x p x u x p x u , x p x v g P l, x p x a g^ ', |x p x^(7^ <J , \xf x^ g v ^ p y g p ^g p ° 



x 



/8 
2 
2 
2 
2 



2 2 2 2 1 \ 

2 10 1 

1 2 1 

2 

2 

1 1 I J 



( 2J[(q T ■ k 1T )(q T ■ k 2T )(q T ■ k 3T )(q T ■ k AT )]/Q^ \ 
-2J[(q T ■ k 1T ){q T • k 2T )(k 3T ■ k 4T )]/Q 2 T 
-2J[(q T ■ k 3T )(q T ■ k 4T )(k lT ■ k 2T )]/Q 2 T 
-J[(q T ■ ki T )(q T ■ k 3T )(k 2T ■ fc 4T ) + 1 <-> 2}/Q 
-J[(q T ■ k 1T )(q T ■ k AT )(k 2T ■ k 3T ) + 1 <-> 2]/Q! 
V2J[(fc 1T -fc 2T )(fc 3 T-fc 4 T)] / 



(A10) 



For the implementation of these rather lengthy formulas we used FORM 25 
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FIGURES 

FIG. 1. Quark and antiquark handbag diagrams for inclusive DIS. 
FIG. 2. The quark and antiquark Born diagrams for the Drell-Yan process. 
FIG. 3. The blob representing the quark correlation function. 
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